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A  ZERO-ONE  INTEGER  PROGRAM  TO  SOLVE  MULTIACTIVITY 
MULTIFACILITY  CAPACITY-CONSTRAINED  ASSIGNMENT 
PROBLEMS 


by 

Krishan  L.  Chhabra 
Richard  M.  Soland 


1 .  Introduction 


Multiactivity  multifacility  assignment  problems  are  important 
because  of  their  practical  applications.  The  objective  of  such  a  problem 
is  to  minimize  the  total  cost  of  the  system,  which  includes  both  variable 
costs  and  fixed  costs.  A  problem  which  has  been  of  recent  interest  takes 
into  account  capacity  constraints  as  well  [Gross,  Pinkus,  and  Soland  (1979)] 
A  solution  algorithm  and  a  computer  program  implementing  it  have  been 
developed  for  this  kind  of  problem,  i.e.,  for  multiactivity  multifacility 
0-1  assignment  problems  with  capacity  constraints. 

The  development  of  the  solution  algorithm  and  computational  test 
results  using  the  computer  program  are  described  in  a  separate  report 
titled  "Solving  a  Multiactivity  Multifacility  Capacity-constrained  0-1 
Assignment  Problem." 


This  document,  on  the  other  hand,  provides  a  description  of  the 
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computer  program  and  instructions  for  its  use.  It  is  a  FORTRAN  program  and  is 
named  ZIPCAP,  an  acronym  for  Zero-one  Integer  Program  to  solve  multi¬ 
activity  multifacility  Capacity-constrained  Assignment  Problems. 

Section  2  of  this  document  provides  a  mathematical  formulation 
of  the  basic  problem  and  examples  of  application  areas.  Section  3 
includes  the  computational  steps  and  a  program  description,  whereas 
Section  4  provides  user  information  including  program  options  and 
input  data  specifications.  Test  problems  and  output  illustrations 
are  covered  in  Section  5. 

2.  Problem  Formulation 


2 . 1  Basic  Problem 


The  basic  problem,  called  problem  (P) 


is  to  find  x. .  and 
i  3 


y^  values  that: 

fMir  imize 

subject  to 


(P) 


m 

n 

P 

E 

i=l 

E 

j=l 

au  “u  +  £  Vk 

(1) 

m 

E 

i=l 

xij 

*  1  j *1 »  • » n 

(2) 

m 

E 

i=l 

n 

E 

j*l 

dijk  xij  -  Vk  kml . . 

(3) 

xij  ’ 

*  0  or  1  for  all  i,  j,  k 

(4) 

V 

where  m  ,  n  ,  and  p  represent  the  number  of  designs,  activities,  and 
facilities  respectively;  and  i  ,  j  ,  k  are  the  corresponding  indices. 
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The  parameters  are  interpreted  as  follows: 


a^  =  variable  cost  of  activity  j  using  design  i  , 
=  fixed  cost  of  facility  k  , 
s^  =  capacity  available  at  facility  k  ,  and 


d 


ijk 


=  capacity  required  at  facility 
when  activity  j  uses  design 


k  for  activity 
i  . 


j 


The  decision  variables  are: 

x„  =1  if  activity  j  uses  design  i  , 

=  0  otherwise. 

=  1  if  facility  k  is  used  , 

■  0  otherwise. 

Simply  stated,  in  problem  (P)  one  seeks  to  minimize  the  sum  of 
variable  and  fixed  costs  subject  to  constraint  set  (2)  that  each  activity 
be  assigned  a  single  design  and  subject  to  constraint  set  (3)  that  the 
capacity  constraints  imposed  by  the  facility  capacities  are  satisfied. 

Problem  (P)  has  mn+p  0-1  variables  and  n+p  constraints. 

2.2  Key  Elements 

The  key  elements  in  problem  (P)  are  activities,  facilities, 
designs,  variable  costs,  and  fixed  costs.  In  general,  these  can  be 
defined  as  follows. 


.  Activity  —  a  type  of  component,  product  or  item  under 
consideration. 
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.  Facility  —  an  installation  or  location  where  activities 
can  be  served  or  serviced. 

.  Design  —  a  meaningful  configuration  of  facilities  along 
with  a  meaningful  strategy  for  using  them.  For  example, 
if  there  are  four  facilities  1,  2,  3  and  4,  a  design 
may  include  facilities  1,  2  and  4,  whereas  another 
design  may  include  facility  2  only. 

.  Variable  Cost  —  the  cost  associated  with  a  particular 
activity  being  assigned  to  a  particular  design. 

.  Fixed  Cost  —  the  cost  associated  with  a  facility 
(such  as  part  of  its  purchase,  operation  and 
maintenance  cost)  that  is  independent  of  the  activities 
served  at  that  facility. 

2.3  Areas  of  Application 

The  computer  program  ZIPCAP  is  designed  to  solve  a  problem  that 
can  be  formulated  as  problem  (P) .  The  formulation  applies  to  existing 
as  well  as  to  proposed  facilities.  In  other  words,  it  is  useful  for  a 
situation  wher<  the  decision  may  be  to  delete  some  of  the  existing 
facilities,  as  well  as  for  a  situation  involving  a  choice  among  proposed 
facilities  (or  locations) . 

Table  1  includes  examples  of  areas  where  this  formulation  can 
be  applied.  Within  each  application  area,  activities  and  facilities 
are  identified.  The  implications  of  designs,  variable  costs  and  fixed 
costs  are  appaient. 

3.  Computational  Procedure  and  Program  Description 

The  procedure  to  solve  problem  (P)  basically  involves  branch 
and  bound  [Ceoffrion  and  Marsten  (1972)],  and  Lagrangean  relaxation 
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[Geoffrion  (1974)].  In  addition  certain  heuristics  are  used  for 
branching,  and  to  exclude  infeasible  assignments. 

Branching  is  done  on  the  x„  variables  by  setting  them  equal 

to  1  or  0  .  The  branching  commences  by  setting  a  variable  equal  to 
1  and  moving  to  the  left-branch  node.  When  backtracking,  the  correspon¬ 
ding  variable  is  set  equal  to  0  and  we  move  to  the  right-branch  node 
(if  the  right-branch  lias  not  yet  been  explored).  Figure  1  illustrates 
a  branching  tree. 


The  x. .  variables  which  have  been  set  equal  to  1  or  0  at 
ij 

any  given  node  .•re  termed  fixed  variables,  and  the  remaining  ones  are 
termed  free  variables. 

3.1  Computational  Steps 

Figure  2  presents  a  simplified  flow  diagram,  showing  the 
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Figure  2 

Simplified  Flow  Diagram  Showing  Computational  Steps 


computational  steps  of  the  computer  program.  Following  is  a  description 
of  these  steps. 


Step  1.  Initialize  by  setting  best  upper  bound  (BUB)  equal  to  a 
large  value  and  the  node  number  to  1. 

Step  2.  Apply  the  capacity  rule  to  exclude  infeasible  assignments 
prior  to  solv  ng  the  relaxed  problem. 


Exclude  a  free  x^  if,  for  any  facility  k  end  activity  j  , 


(dijk  -  mjn  dijk) 

> 

(s* 

-  •}”  dijk)  . 

If  an  : 

\i..  is  already  fixed 
ij 

as  1 

,  the  corresponding  d  k 

replaces  min 

j 

d. . 
ijk 

If  Zj 

min  dijk  > 

\  • 

then 

backtrack,  i.e.,  go  to  step  12. 

Step  3. 

Solve  the  relaxed  problem 

(PR)  ,  i.e., 

f 

Minimize 

m 

Z 

i=l 

n 

Z 

3-1 

Cij  Xij 

(PR)  - 

subject  to 

m 

Z 

i=l 

XU 

=  1  j=l  ,  •  •  •  ,  n  , 

Xij 

=  0  or  1  for  all  i,  j  * 

where  i 

cij '  + 

P 

Z 

k-1 

bk  (] 

1  "  4k  )  (dljk/ok). 

6k-  i  «  i 

5 

XiJ 

eik  >  0  • 

such  that 
x^  is  fixed 


0  otherwise,  and 
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=  1  if  design  i  uses  facility  k  , 

=  0  otherwise. 

Step  4.  Lower  bound  (LOWB)  is  the  sum  of  the  optimal  solution 

value  of  problem  (PR)  and  the  fixed  cost  FC.  FC  is  the  cost  of  the 

facilities  forced  into  the  solution  because  of  the  x..  variables 

ij 

fixed  as  1: 


FC 


P 

E 

k=l 


Vk 


Step  5.  Compare  lower  bound  with  best  upper  bound.  If  LOWB 
>  BUB  ,  go  to  step  12. 

Step  6.  Check  if  the  solution  of  problem  (PR)  obtained  in 
step  3  satisfies  the  capacity  constraints  of  problem  (P)  ,  i.e.. 


m 

E 

i=l 


n 

E  d 
j=l 


X  < 

ijk  xij  - 


skYk  k=l,  •••*  P» 


where  y,  =  1  if  E.  E,  d.,,  x.  .  >  0  , 

k  i  j  ijk  ij 


=  0  otherwise 


If  the  capacity  constraints  are  not  satisfied,  go  to  step  10. 


Step  7.  Compute  upper  bound  (UPB)  ,  i.e., 

m  n  p 

UPB  ■  E  E  a  x,  +  E 

i-1  j=l  J  3  k-1 


Vk 
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Step  8.  Compare  upper  bound  with  best  upper  bound.  If  UPB  >_  BUB, 
go  to  step  10. 

Step  9.  Retain  UPB  as  BUB  ,  i.e.,  set  BUB  =  UPB  and  note 
the  corresponding  solution  comprising  x^  and  variables. 

If  I.OWB  equals  BUB  (same  as  UPB  in  this  case),  go  to  step  12. 


Step  10  If  an  variable  has  been  fixed  as  1  for  each  of  the 

n  activities,  then  go  to  step  12. 


Step  11.  Select  an  x„  variable  for  branching  to  the  left  node. 
The  choice  of  the  branching  variable  depends  on  the  c^  values  and  is 
such  that  the  corresponding  x„  ,  if  perturbed,  has  the  maximum  impact 
on  the  optimal  value  of  problem  (PR).  That  is,  for  each  j  ,  obtain 
Dj  =  c^  -  c^  where  c^  is  the  second  smallest  permissible  and 
is  the  smallest  permissible  c  value  in  column  j  . 


A  c^  value  is  considered  permissible  if  it  does  not  correspond 

to  an  activity  j  having  x^  variable  fixed  as  1  ,  or  to  an  x^ 

variable  fixed  as  0  .  The  x..  variables  can  be  fixed  as  1  or  0 

ij 

because  of  the  branching  rule,  backtracking,  and/or  use  of  the  capacity 
rule.  Further,  may  be  fixed  as  0  because  of  the  bounding  rule. 


“  (c«  '  c«)  >  (' 


BUB  -  LOWBj  . 

The  x. .  variable  selected  has  the  maximum  value  of  D.  and  has  c 


ij 


c, , 


lj  * 


Set  the  selected  x^  variable  to  1  .  Increase  the  node  number 
by  one,  and  go  to  step  2. 


Step  12.  Backtracking  continues  until  a  node  is  reached  for  which 
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the  right  branch  has  not  been  explored.  Set  the  corresponding  x 
variable  to  0  .  Increase  the  node  number  by  one,  and  go  to  step 


ij 


2. 


If  no  such  node  is  found,  it  implies  that  the  root  node  has  been 
reached,  and  the  procedure  terminates. 

The  current  BUB  and  the  solution  obtained  in  step  9  provide 
the  optimal  s(  lution  and  its  value.  If  BUB  happens  to  be  the  initial 
value,  the  problem  does  not  have  a  feasible  solution,  i.e.,  the  capacity 
constraints  cannot  be  satisfied. 

3.2  Computer  Program 

The  computer  program  ZIPCAP  is  written  in  FORTRAN  IV.  It  is 
based  on  the  flow  diagram  and  the  computational  steps  described  in  the 
preceding  paragraphs. 

Appendix  A  provides  a  listing  of  the  program.  Extensive  use 
of  comment  cards  provides  details  within  each  computational  step. 

A  dictionary  of  the  symbolic  names  used  in  the  program  listing 
is  provided  in  Appendix  B. 

The  program  has  been  developed  and  tested  on  an  IBM  3031 

is 

at  The  George  Washington  University.  The  only  internally  supplied 
subroutine  used  by  the  program  is  TIMET  from  the  PI)  Library.  This 
subroutine  provides  the  elapsed  time  for  program  execution  and  excludes 
the  time  to  read  and  write  the  input  and  output. 


* 

IBM  3031  CPU  vith  3  million  Bytes  Real  Memory.  Operating  System 
0S/VS1  Release  6.7.  Various  compilers  including  FORTRAN  levels 
G  1  to  6  and  H  1  to  6. 
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The  program  comprising  about  480  lines  is  currently  dimensioned 
for  a  problem  size  of  35  designs  (m)  ,  35  activities  (n)  and  30 

facilities  (p)  .  The  user  capacity  required  to  execute  this  program 
is  346  K  bytes.  The  functional  relaw ionship  between  the  dimensions 
of  the  arrays  and  the  storage  requirement  is  given  by 


f  (m  ,  n 


p)  -  4 


[< 


(p+4)mn  +  (m+5)p  v 


bytes . 


Thus  the  program  size  to  execute  a  problem  has  two  components: 
one,  due  to  the  program  itself,  comprises  173  K  bytes,  and  the  other, 
dependent  on  the  dimensions  of  the  arrays,  is  given  by  the  above 
functional  relationship. 


4.  User  Information 


ZIPCAP  is  designed  to  solve  a  multiactivity  multifacility 
capacity-constrained  0-1  assignment  problem  having  formulation  (P)  . 
Such  a  problem  is  referred  to  as  "capacitated"  because  the  optimal 
solution  murt  satisfy  the  capacity  constraints.  However,  the  program 
can  be  used  in  a  situation  where  the  capacity  constraints  are  not 
applicable,  in  other  words,  for  an  "uncapacitated"  problem.  The  program, 
in  this  case,  will  generate  the  necessary  informat_on  to  conform  to 
formulation  (P)  . 

Figure  3  presents  a  schematic  diagram  of  the  deck  structure 
for  using  ZfPCAP.  The  cards  include  both  the  job  control  cards  and  the 
input  data  c; rds  required  to  solve  a  problem. 

Following  is  a  detailed  description  of  these  cards  (in  the  same 
order  as  presented  in  Figure  3);  it  provides  the  necessary  coding 
information. 

(1)  Job  and  JCL  Cards:  Figures  4A  and  4B  show  two 
alternatives  for  the  use  of  JCL  cards. 
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(9)  End  Card 


ZIPCAP  Data  Structure 


Job  and  JCL  Cards  —  Alternative 


FORTRAN  STATEMENT 


Job  and  JCL  Cards  —  Alternative  II 
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Altei native  I  uses  the  source  program  for  the  compilation  and 
execution  procedure.  In  Alternative  II,  the  program  is  first  compiled 
and  stored,  and  thereafter,  for  a  given  problem,  onlv  the  execution  is 
necessary. 

The  program  uses  the  TIMET  subroutine  from  the  PL1  Library 
in  order  tu  record  the  elapsed  time.  This  time  includes  only  the 
execution  time  of  the  program  and  excludes  the  time  needed  to  read  the 
input  and  to  write  the  output.  The  JCL  cards  needed  for  this  purpose 
are  identified  in  Figures  4A  and  4B  .  If  this  subroutine  is  not 
available,  the  corresponding  JCL  cards  are  not  requi  ed.  In  that  case 
the  appropriate  source  program  cards  should  also  be  removed.  (These 
cards  are  identified  in  the  program  listing  in  Appendix  A.) 

(2)  Options  Card 

The  options  card  is  the  first  input  data  card  and  is  described 
in  Table  2. 


(3)  Parameter  Card 

Table  3  describes  the  parameter  card. 

Table  3 

THE  PARAMETER  CARD 


Columns 

Format 

Name 

Definition 

1-5 

15 

M 

Number  of  designs  (m) 

6-10 

15 

N 

Number  of  activities  (n) 

ll-j.5 

15 

P 

Number  of  facilities(p) 
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Even  with  ISTEP  =  0,  the  output 
includes  the  total  number  of  no 
explored . 


Columns  Format  Name  Value  and  Definition  Remarks 
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(4)  to  (8)  Other  Input  Data  Cards 

Table  4  describes  the  cards  for  the  variable  costs  a^.  , 
the  fixed  costs  b^  ,  the  capacity  availabilities  s^  ,  and  the 
capacity  usage  values  d^^  •  However,  if  solving  an  uncapacitated 

problem,  i.e.,  where  capacity  constraints  are  not  active,  the  facility- 
design  values  e^  are  required  as  input  data  instead  of  the  s^  and 
d^jk  values. 

The  input  data  values  for  a,.  ,  b.  ,  s,  .  and  dJM  are 

J-3  k  k  ijk 

required  ii.  integer  format.  Any  real-valued  data  can  be  converted  into 
integer  values  by  scaling,  and  the  optimal  solution  reconverted  to  real 
valyes  later.  For  example,  real-valued  data  of  215.35  ,  116.50  ,  .... 

180.61  can  be  treated  as  21535  ,  11650 .  18061  ;  and  an  optimal 

value  of  584230  obtained  from  the  program  would  corr< spond  to  5842.30. 

(9)  End  Card 

This  is  the  usual  end  card  following  input  data  cards.  It 
contains  '//'  in  the  first  two  columns. 

5.  Test  Problems  and  Output  Illustrations 

5.1  Test  Problem  1 

5.1.1  Problem  formulation  (stated  for  ease  of 
reference) 


The  problem  is  to  find  x^  and  y^  variables  that 
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OTHER  INPUT  DATA  CARDS 


IUNCAP=1,  this  informati 
generated  by  the  program. 


Table  4 
(Continued) 


uncapacitated  problem, 
i.e.,  for  IUNCAP=1 . 


Minimize 


subject  to 


where  m  =  3 


n 

m 

P 

Z 

3“1 

Z 

i=l 

a  .  .  x. . 
>■3  ij 

+ 

1  bkyk 
k=l  *  K 

m 

Z 

x. . 

=  1  J 

l-i. 

•  •  i  n 

i=l 

ij 

n 

m 

I 

z 

d  ,  . ,  x .  , 

< 

s,  y,  k=l , . 

j=l 

i=l 

ijk  ij 

k^k  ’ 

Xij 

,  yk=  0 

or 

1, 

n  =  4  ,  p  *  5. 


The 


a.,  values  are: 
ij 

j 

1  2  3  4 


bx  -  1750  b2  =  2000  b3  =  1750  b 4  *  1350 

s,  =  400  s„  =  400  s_  =  1000  s.  =  400 

1  2  3  4 


The  d  .  . . 


values  are  as  follows 
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5.1.3  Annotated  Output:  The  computer  output  for  this  test 
problem  is  presented  in  Appendix  C.  Most  segments  of  this  output  are  self- 
explanatory.  The  order  of  these  segments  is  lettered  in  circles  as  follows 


(a)  Options  selected  —  information  given  on 
options  card 

(b)  m,  n,  and  p  —  information  from  parameter  card. 
This  segment  is  printed  if  the  option  IINPT=1 


(c)  Other  input  data,  i.e.,  a^  ,  b^  ,  s^  and 

d  ^  for  a  capacitated  problem;  and  as 

generated  by  the  program.  For  an  ’incapacitated 
problem,  input  data  a^  ,  b^  ,  and  e^  ; 
and  and  d^^  as  generated  by  the  program. 

This  segment  is  printed  if  the  option  IINPT=1 

(d)  Intermediate  steps'  display,  depending  on  the 
value  of  the  option  ISTEP.  No  display  for 
ISTEP=0,  summary  for  ISTEP=1,  and  detailed 
steps  for  ISTEP=2 

(e)  Elapsed  time  in  seconds  to  execute  the  program, 
excluding  the  time  to  read  the  input  data  and 
to  write  the  input  and  output.  It  includes 
the  time  to  print  intermediate  steps  if 
option  ISTEP=1  or  2 

(f)  Total  number  of  nodes  explored  for  the  problem 
solved 

(g)  Optimal  solution  showing  the  x^  variables 

with  value  1,  y^  variables  with  value  1  or  0, 

and  the  value  of  the  objective  func“ion.  If  the 
problem  does  not  have  a  feasible  solution,  this 
segment  prints  "Problem  does  not  have  a  feasible 
solution  because  the  capacity  constraints  cannot 
be  satisfied." 


The  node  and  bounds  information  at  a  specified  time  is  displayed 
only  if  option  FT  is  assigned  a  particular  value. 
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5.2  Test  Problem  2 


This  is  a  capacitated  problem  having  5  designs,  4  activities  and 
8  facilities. 


The  problem  formulation  and  coded  input  are  similar  to  that  of 
test  problem  r. 

The  annotated  output  (which  includes  the  listing  of  input  data) 
is  presented  in  Appendix  D.  The  annotations  are  the  same  as  those 
already  described  for  test  problem  1.  Note  that  the  options  selected 
include  ISTEP=0,  hence  no  intermediate  steps;  and  EPS=0.002,  implying 
that  the  solution  may  be  suboptimal  within  +  0.2%  . 

5.3  Test  Problem  3 


Although  ZIPCAP  is  designed  for  capacitated  problems,  it  may 
be  used  to  solve  uncapacitated  problems  as  a  special  case.  For 
demonstration,  this  test  problem  is  uncapacitated,  having  10  designs, 
8  activities,  and  8  facilities.  The  input  data  include  the  options 
card,  the  parameter  card,  input  data  values  for  a^  ,  b^  ,  and 


e 


ik  ' 


The  program  generates  the  d  ^  and  s^  values. 


Appendix  E  shows  the  annotated  output  for  this  problem.  Note 
that  the  options  selected  include  IINPT=1,  hence  the  data  listing; 
ISTEP=0,  therefore  no  intermediate  steps;  ICAPR=0  and  IUNCAP=1  being 
an  uncapacitated  problem.  Note  that  option  ET  has  been  specified  a 
value  of  3C.0  .  Thus,  if  the  optimal  solution  was  not  reached  within 
30  seconds,  the  node  and  bounds  information  at  time  ET  would  have  been 
displayed . 
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MAIN 


DATE  »  00206  14/30/07 


0001 


c 

L 

c 
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0002 

0003 

0004 

0009 
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0007 


0000 
C  0004 

.  0010 
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,  0012 

C  0013 


0014 


0019 

0016 

0017 

0010 

0014 

0020 

0021 

0022 


c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 
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c 

c 

c 
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Z1PCAP ,  A  ZERO-ONE  INTEGER  PROGRAM  IS  DESIGNED 
TO  SOLVE  MULTI ACT IVITV  MULTIFAC1L1TV  CAPACITY- 

CONSTRAINED  problems  having  variable  and  fixeo 
COSTS.  IT  ALSO  SOLVES  UNCaPaCITaTED  PROBLEMS  AS  A 
SPECIAL  CASE 

INTEGER  D(35,35,30),  AI39.39I,  CX(35,35),  E(35,30>, 

1  6130) ,  BS0LXI35) «  BSOLVOO),  FLB ( 30) ,FIX(35 > . FIX1 135 ) . 

2  FUB ( 30 ) «  S  (30)  ,  S0LX(35),  STXI1225) 

REAL  MINC  1 35 )  t  NMINCOS) 

DIMENSION  C(35,35>,  DIFBR<35>,  XT2I35),  M1NDI35) 

INTEGER  6R0.  BR1»  PC.  FCUB,  P 

REAL  LOWBt  MAXD1F.  M1NSC 

•»»,******»*««0PT1UNS  available:  IINPT,  ICAPP,  ISTEP,  I UNCAP t EPS 
IINPT-1  IF  INPUT  LISTING  DESIRED;  0  OTHERWISE 
IC APR- 1  IF  CAPACITY  RULE  TO  BE  USED;  0  OTHERWISE 
ISTEP-0  IF  LISTING  OF  INTERMEDIATE  STEPS 
NOT  DESIRED.  ISTEP-1  IF  SUMMARY  OF  BRANCH  t 
BOUNO  NODES  DESIRED.  ISTEP-2  IF  DETAILED 
LISTING  OF  INTERMEDIATE  STEPS  DESIRED. 

1UNCAP-1  IF  SOLVING  AN  UNCAPACITATED  PROBLEM, 

0  OTHERWISE. 

EPS-  A  FRACTIONAL  VALUE  IF  SuBOPTINAL 
SOLUTION  0ES1RE0,  E.G.,  EPSILON  AS  0.005 
IMPLIES  SOLUTION  TO  BE  WITHIN  -0.5  PERCENT 
OF  THE  OPTIMAL  SOLUTION.  EPS-0. 0  IF  OPTIMAL 
SOLUliOF  DESIRED. 

ET-  ELAPSED  TIME  IN  SECONDS,  IF  SPECIFIED.  AT 
WHICH  THE  NODE  AND  BOUNDS  RELATED  INFORMATION 
IS  PRINTEO.  THIS  IS  USEFUL  IN  A  SITUATION  IF 
1  STEP-0  AND  THE  PROGRAM  TERMINATES  BEFORE 
REACHING  The  FINAL  SOLUTION. 

•••••••••••t—RtAO  input  data*-—— •——  ———*  — 

READ  10,  IINPT,  1CAPR,  I  STEP,  IUNCAP,  EPS,  ET 
10  FORMAT  (Nil,  F6.5,  F10.3) 

M-  NUMBER  OF  OESIGNS 
N-  NUMBER  OF  ACTIVITIES 
P-  NUMBER  OF  FACILITIES 

RE-0  20,M,N»P 
20  FORMAT  (315) 

A(I.J):  VARIABLE  COST  MATRIX 
READ  30,  ((A(I,J),  I- 1 ,M ) , J-l ,N ) 

30  FORMAT  (0110) 

BIX):  FIXED  COST  VECTOR 
REAO  30,  (0(K)»X-1,P) 

IF  I IUNCAP. EO. 1)  GU  TO  40 

S(K) I  CAPACITY  LIMIT  VECTOR)  REGUlREO  ONLY 
IF  IuNCAP-0 
REAO  30,  (S(K),K-1,P> 

D(l,J,K)l  CAPACITY  USAGE  MATRIX;  REQUIRED 
ONLY  IF  IUNCAP-0 

DO  32  K-l.P 

READ  30, ((0(1, J, A),  1- 1 ,M I , J-l »N ) 

32  COMlNuE 
DO  37  X-l.P 
00  37  I • l ,M 

IF  (0(I,1,K).EQ.0I  GO  TO  35 

E(I,AI*1 

GO  TO  37 


00000010 

00000015 

00000020 

00000030 

00000040 

00000050 

00000060 

00000070 

ooooooeo 

00000040 

00000110 

00000120 

00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000140 

00000200 

00000210 

00000220 

00000230 

00000240 

00000250 

00000253 

00000256 

00000260 

00000263 

00000266 

00000270 

00000273 

00000200 

00000240 

00000300 

00000310 

00000320 

00000330 

00000360 

00000350 

00000360 

00000370 

00000300 

00000340 

00000600 

00000610 

00000420 

00000430 

00000440 

00000450 

00000660 

00000470 

00000600 

00000440 

00000900 

00000910 

00000920 

00000530 


30  - 
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0023 

35  EII,KI*0 

00000540 

0024 

37  CONTINUE 

00000550 

0025 

GO  TO  90 

00000560 

C 

E(I,K|I  DESIGN-FACILITY  MATRIX!  REQUIRED  ONLY 

00000570 

C 

IP  I UNCAP*  1 

00000580 

0026 

40  REAO  45,(U(I,X),1*1,M>,X*1,PI 

00000590 

0027 

45  FORMAT  (80111 

00000600 

0028 

00  80  X*1,P 

00000610 

0029 

S(K»«N 

00000620 

0030 

00  75  I-1,N 

00000630 

0031 

IF  UdtM.EQ.il  GO  TO  65 

00000640 

0032 

00  60  J>ltN 

00000650 

0033 

D(I,J,Kl-0 

00000660 

0034 

60  CONTINUE 

00000670 

0035 

GO  TO  75 

00000680 

0036 

65  00  70  J*1,N 

00000690 

0037 

01 1  »J»M*1 

00000700 

0038 

70  CONTINUE 

00000710 

0039 

75  CONTINUE 

00000720 

0040 

80  CONTINUE 

00000730 

c 

••••••••••••••PAINT  INPUT  OAT A ••••••••*•••••••••••••• •••••*••• 

00000740 

0041 

90  PRINT  95,  1INPT,  ICAPR,  l&TEP,  IUNCAP,  EPS,  ET 

00000750 

0042 

95  FORMAT  1  •!• , •  OPTIONS  SELECTEO  I  1INPT*‘,I1, 

00000760 

1  •  ICAPR* *,11,  •  ISTEP** ,11,  •  IUNCAP**, 11, 

00000770 

2  •  EPS** ,F8. 5,  •  ET**,  F10.3///I 

00000780 

0043 

IF  (UNPT.EU.OI  GO  TO  168 

00000790 

0044 

PRINT  1 00  .  M  ,N  ,P 

00000800 

0045 

100  FORMAT  (*0*,  T55,  ‘INPUT  DATA*, /IX,  T55,  . — - -  ——»,////// 1 X  , 

00000810 

1T41,  •NUMBER  OF  DESIGNS  (Ml**,  4X , 14//1X,T41 , 

00000820 

2 ’NUMBER  OF  ACTIVITIES  (Nl**,  1X.14//1X,  T41, 

00000830 

3 ‘number  of  facilities  ipi**,1x,  in/z/i 

00000840 

0046 

►  PINT  105 

00000850 

004.7 

105  FORMAT  (  4 X ,  ‘VARIABLE  COST  MATRIX  All,J|‘,/4X, 

00000860 

00000870 

0048 

00  110  1*1, M 

00000880 

0049 

110  PRINT  115,  I,  ( A ( I ,JI, J- 1 , N 1 

00000890 

0050 

115  FORMAT  (*0‘,  T6,  •  !••»  13,  4X,81 13,  41/,  1*-X,BU3II 

00000900 

0051 

PRINT  120 

00000910 

0052 

120  F0RMAT(‘0‘,//4X,‘FIXE0  COST  VECTOR  B(Rl‘,/4X, 

00000920 

00000930 

0053 

PRINT  122,  (81X1 ,X*1 ,P 1 

00000940 

0054 

122  FORMAT  (‘0‘,  T 15 ,  8113,  3(/,  14X, 811311 

00000950 

0055 

PRINT  125 

00000960 

0056 

125  FORMAT (• O', //4X, ‘CAPACITY  LIMIT  VECTOR  S(Rl*,/4X, 

00000970 

00000980 

0057 

PRINT  128,  ISIM,K-1,PI 

00000990 

0058 

128  FORMAT  l‘0‘,  T15,  8113,  3(/,  14X.8U3II 

00001000 

0059 

PRINT  130 

00001010 

0060 

130  FORMAT! *0*  »//4X ,  ‘CAPACITY  USAGE  MATRIX  0(1, J, XI*  »/4X» 

00001020 

1’-— - - - .-.  —  ■»,/  | 

00001030 

0061 

00  150  X*1,P 

00001040 

0062 

PRINT  135,X 

00001050 

0063 

135  FORMAT  ( ‘0‘,//5x,‘X*‘,I3/l 

00001060 

0064 

00  145  1*1 ,M 

00001070 

0065 

PRINT  140,1, (0(1, J, XI,  J— 1 ,N 1 

00001080 

0066 

140  FORMAT  (‘0‘,  T 6,  •!*••  13,  4X,SI 13,  4(/,  14X,8I13II 

00001090 

0067 

145  CONTINUE 

00001100 

0068 

150  CONTINUE 

00001110 

31 
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PRINT  155 

i  FORMATt'0',//4X, 'DESIGN-FACILITY  MATRIX  EII,K)',/4X, 

l—  — . —  ',/) 

00  160  I-l.M 

PRINT  15b,  I,  (Efl.K ),K>1,P) 

I  FORMAT  1*0*,  T6,  •!••»  13,  4X,SI13,  31/,  14X.8I13)) 

)  CONTINUE 

)  IF  I1STEP.EQ.0)  GO  TO  190 
IF  I ISTEP.EQ. 1 )  GO  TO  175 
PRINT  170 

1  FORMAT  (  •0*,///55X, • DETAILED  LISTING  OF  STEPS',/) 

GO  TO  190 
)  PRINT  180 

)  FORMAT  I  '0 1 , ///55X •' SUMMARY  OF  STEPS',/) 

)  8UB«9999999. 

BUBS-  BUB/  ll.O'EPS) 

NSX-0 
NOO-l 
18M00-1 
INPT-0 
INSET-0 
00  205  J«1 ,N 
FIX  I J )— 0 
KT2IJ1-0 
00’ 205  1-1, M 
CX«I,J)-0 
CONTINUE 
LQ1-0 
LQ2-0 
LR2-0 

CALL  TIMETIITO) 

1FINSX.EQ.0)  00  TO  283 

CXIl.J)  CONTAINS  fixed  AND  FREE  XII, J)  VARIABLES. 
STXIINS)  CONTAINS  FIXED  XII, J)  VARIABLES. 

CXII,J)  AND  STXIINS)  ARE  UPDATED  BY  THE  CAPACITY 
RULE,  THE  BOUNDING  RULE,  ANO  THE  RULE  FOR 
BRANCHING  AND  B AC XT RACK I NO. 

IN  CXI  I , J )  A  FIXED  VARIABLE  IS  RECOROED  AS  1  OR 
2,  ANO  A  FREE  VARIABLE  AS  O. 

A  VALUE  OF  1  IMPLIES  THAT  THAT  PARTICULAR  VARIABLE 
IS  FIXED,  AND  FIX(J)  IS  SET  EQUAL  TO  1  IMPLYING 
THAT  COLUMN  J  HAS  A  FIXED  VARIABLE  OF  VALUE  1. 

A  VALUE  OF  2  IMPLIES  THAT  THAT  PARTICULAR  VARIABLE 
SHOULO  NOT  BE  CONSIDERED  FOR  CURRENT  COMPUTATIONS. 
AN  Xll,jl  RECOROEO  IN  CXti.J)  AS  I  DUE  TO  THE 
BRANCHING  RULE  IS  RECORDED  IN  STXIINS)  AS  X«100«J. 
AN  XII, J)  RECORDED  IN  CX(1,J)  AS  1  DUE  TO  THE 
CAPACITY  RULE  OR  THE  BOUNDING  RULE  IS  RECORDED  IN 
STXIINSI  AS  I  X* 1 00a J ) *  1000000 . 

AN  XII, J)  RECORDED  IN  CXll,U)  AS  2  IS  RECORDED  IN 
STXtINSi  AS  •IX-100«J >•1000000. 

IF  IISTEP.EO.OI  GO  TO  225 
PRINT  220, NOD 

FORMAT  I'O',//  6X«  * NOOE  NUMBER', 15/) 

••••*•  ••••P-9-UP0ATE  CXtI,J|  FOR  BRO*--—--***— 

BRO  IS  THE  RIGHT  BRANCHING  VARIABLE 

LX-BRO 


00001120 

00001130 

00001140 

00001150 

00001160 

00001170 

00001160 

00001190 

00001200 

00001210 

0000122C 

00001230 

00001240 

00001250 

00001260 

00001270 

00001280 

00001290 

00001310 

00001315 

00001320* 

00001330' 

00001390 

00001400 

00001410 

00001420 

00001430 

00001433 

00001436 

00001440 

00001443 

00001445* 

00001447 

00001450 

00001460 

00001460 

00001490 

00001500 

00001505 

00001510 

00001515 

00001520 

00001525 

00001530 

00001535 

00001540 

00001545 

00001550 

00001555 

00001560 

00001565 

00001570 

00001580 

00001590 

00001600 

00001610 

00001615 

00001620 


*Relevant  to  the  TIMET  Subroutine 
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0104 

IX*LX/100 

00001630 

0105 

JX*LX-IX*100 

00001640 

0106 

CXIIX, JXI*2 

00001650 

0107 

KT2I JXI-KT2I  JX)*1 

00001660 

0108 

FIX! JXI*0 

00001720 

0109 

LQ1=LQ1“1 

00001725 

0110 

IF  IKT21 JX).LT.(»»1I 1  GO  TO  270 

00001730 

0111 

DO  255  2  *1  *N 

00001740 

■_ 

0112 

IF  (CXII.JXl.EQ.2l  GO  TO  255 

00001750 

0113 

CXIli JX1*1 

00001760 

0114 

NSX*NSX+1 

00001763 

0115 

STXINSXI-  CI*100*JX >♦1000000 

00001766 

0116 

FIX1 JX>*1 

00001770 

0117 

L01>L0I*1 

00001780 

0118 

FI  XI ( JX 1  =  I 

00001790 

0119 

GO  TO  270 

00001800 

0120 

255 

CONTINUE 

00001810 

0121 

270 

LQ2*0 

00001820 

0122 

LR2*0 

00001825 

0123 

GO  TO  283 

00001830 

V 

0124 

272 

IF  (ISTEP.EQ.OI  GO  TO  276 

00001840 

0125 

PRINT  220. NOD 

00001850 

C 

****«» »*»**4**UP0ATE  CXII.Jl  FOR  BR1***‘ 

00001853 

C 

BR1  IS  THE  LEFT  BRANCHING 

VARIABLE 

00001856 

0126 

276 

LQ2«0 

00001860 

0127 

LR2*0 

00001866 

v/ 

0128 

LX=8R1 

00001870 

0129 

1X-LX/100 

00001875 

0130 

JX=LX-IX*100 

00001880 

0131 

CXC IX. JX 1*1 

00001885 

0132 

FIX! JXI*1 

00001890 

0133 

LQl*LQld 

00001892 

0134 

DO  279  I*1,M 

00001895 

0135 

IF  IIX.EQ.il  GO  TO  281 

00001697 

0136 

279 

CONTINUE 

00001900 

0137 

281 

FlXli JX ) *IX 

00001902 

0138 

283 

IF  (ISTEP.NE.2)  GO  TO  300 

00001905 

( 

0139 

28  5 

DO  295  1*1, M 

00001910 

0140 

PRINT  290.  1.ICXI I. Jl ,J*1,NI 

00001920 

0141 

290 

FORMAT  I  /5X  ,  *CX (I,J)',4X,,I**,I3,2X, 

2014/23X,  20141 

00001930 

0142 

295 

CONTINUE 

00001940 

0143 

PRINT  297,tFIXIJI,J*l,NI 

00001950 

0144 

297 

FORMAT  l/5X,,FIXIJ),,12X,  20I4/23X,  20141 

00001960 

C 

c 

00001970 

c 

AND  UPDATE  CXII.Jl  AND  STXIINS). 

00001980 

0145 

300 

IF  (IUNCAP.EO.il  GO  TO  2300 

00002010 

0146 

305 

IF  I ICAPR.EU.OI  GO  TO  2300 

00002020 

0147 

310 

DO  2000  *»1,P 

00002030 

c 

FINO  THE  SUM  OF  MINIMUM  DII.J.K)  OVER 

EACH  J  FOR  A 

00002040 

c 

GIVEN  X.  I.E.,  HINSD*  SUM 

OF  MINOIJI 

00002050 

0148 

MINSD*0 

00002060 

0149 

DO  900  J*1,N 

00002070 

0150 

IFt  F1X( J I . EU. 0)  GO  TO  350 

00002080 

c 

IF  FIX(JI*1,  SET  M I NO  I J 1 1 

"OII.J.K)  FOR 

CXII.Jl*! 

00002090 

c 

AND  MOVE  10  NEXT  COLUMN  J 

00002 100 

0151 

INO t*F IX 1 ( J 1 

00002110 

0152 

MI  NO! J )*D(lN01,J«Kl 

00002120 

0153 

GO  TO  800 

00002130 
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0154 

350 

LK*0 

00002160 

0155 

1*1 

00002170 

S.- 

0156 

HINDI J1=DI l.J.Kl 

00002180 

C 

SKIP  OII.J.KI  WHEN  CXI  I .J 1 *2  6  HOVE  TO  NEXT  ROM  I 

00002190 

0157 

400 

IFtCXII.JI.EQ.21  GO  TO  600 

00002200 

J 

0158 

500 

IF  101 I.J.K1.LT.MIN0I Jll  HINDI JI*DII.J.KI 

00002210 

0159 

GO  TO  700 

00002220 

0160 

600 

LK*LK*l 

00002230 

'  w' 

0161 

IFII.GT.LM  GO  TO  700 

00002240 

0162 

I*I*l 

00002250 

0163 

HINDI J 1=01 1  .J  .K  I 

00002261. 

vx1 

0164 

GO  TO  750 

00002270 

0165 

700 

I  *  I  ♦  1 

00002280 

0166 

750 

IFI1.LE.H1  GO  TO  400 

00002290 

o 

0167 

800 

MINSO=HINSO*MINDI  Jl 

00002300 

0168 

900 

CONTINUE 

00002310 

0169 

910 

i  r 

1 ISTEP.NE .2  >  GO  TO  980 

00002320 

u 

0170 

PRINT  950.  K.  H1NS0 , 1  HI  NO IJ).J*1.NI 

00002330 

0171 

950 

FORM ATI *0*. *K,  MI NSO . 1  HI NOI J » , J=1 ,N 1 • , 101 1 0 1 

00002340 

C 

F1N0  BALANCE  AVAILABLE  CAPACITY  1BAL0  FOR  A  GIVEN  K 

00002350 

lw 

C 

IF  IBALO  IS  NEGATIVE.  THEN  oACKTRACK. 

00002360 

0172 

980 

IBAL0*SIKI-M1NS0 

00002380 

0173 

1000 

IF 

IIBALD.LT.01  GO  TO  6200 

00002390 

L 

0174 

00 

1500  J«1,N 

00002400 

C 

SKIP  COLUMN  J  IF  FIXI  Jl*l 

00002410 

0175 

IF  IFIXUI.EQ.il  GO  TO  1500 

00002420 

{ ’ 

w 

0176 

00 

1300  1*1. H 

00002430 

c 

SKIP  ROM  I  IF  CXII.JI-2 

00002440 

0177 

1100 

IFICXII.J1.EQ.2I  GO  TO  1300 

00002450 

‘y-*' 

c 

COMPUTE  DIFFERENCE  BETWEEN  DII.J.Kl  AND  HINDIJ1. 

00002470 

c 

IF  IT  IS  MORE  THAN  AVAlLABE  BALANCE.  SET  CXII»JI*2 

00002480 

0178 

1200 

101 F0*0 1 I . J .K 1”HIN0 1 J 1 

00002490 

0179 

IF  ttIDIFD-IBAL01.LE.0l  GO  TO  1300 

00002510 

0180 

CXI I.J)*2 

00002520 

0181 

NSX*NSX+1 

00002523 

0182 

STXINSX1*-! 1*100* Jl-1000000 

00002526 

c 

LQ2  COUNTS  THE  NUMBER  OF  CXII.JI  VALUES  SET  EQUAL 

00002530 

c 

TO  2  IN  A  CYCLE 

00002540 

L 

0183 

LQ2*LQ2*1 

00002550 

c 

KT2IJ1  KEEPS  AN  ACCOUNT  OF  CXII.Jl  VALUES  SET  EQUAL 

00002560 

C 

TO  2  FOR  COLUMN  J 

000025.0 

t 

0184 

KT2 1 J  l“KT2 1 J 1 ♦! 

00002580 

c 

FOR  COLUMN  Jt  IF  ALL  BUT  ONE  CXII.JI  VALUES  ARE 

00002590 

e 

c 

EQUAL  TO  2,  SET  THAT  CX|I.JI*1  6  SET  FIX|J|«1 

00002600 

0185 

1FIKT2IJ1.LT. IH-1 II  GO  TO  1300 

00002610 

0186 

00 

1250  LR*1.M 

00002620 

i„ 

0187 

1FICXILM.JI.EQ.2I  GO  TO  1250 

00002630 

0188 

CXI LR • J 1*1 

00002640 

0189 

NSX*NSX*1 

00002643 

0190 

STXINSX1*  ILR*1 00*J l* 1000000 

00002646 

0191 

F1XIJ1-1 

00002650 

c 

LQ1  KEEPS  AN  ACCOUNT  OF  COLUMNS  FOR  WHICH  FIXIJ|*1 

00002655 

0192 

L01*L01*1 

00002660 

c 

FIXI 1 J 1  SPECIFIES  INDEX  1  FOR  WHICH  FIXIJl*! 

00002662 

0193 

FIXIIJI-LR 

00002665 

0194 

GO  TO  1500 

00002670 

.  ■ 

0195 

1250 

CONTINUE 

00002680 

0196 

1300 

CONTINUE 

00002690 
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0197 

1500 

CONTINUE 

00002700 

0198 

1800 

IF  (ISTEP.NE.2)  GO  TO  2000 

00002710 

0199 

PRINT  1900,  A,  L02 ,  L01 

00002720 

0200 

1900 

FORMAT  1  *0*  ,  *K**  ,13,  •  LG2*',I3,  •  LQ1*‘,  131 

00002730 

0201 

DO  1930  1*1, M 

00002790 

(w. 

0202 

PRINT  290,  1,<CX(I,J),J*1,N> 

00002750 

0203 

1930 

CONTINUE 

00002770 

c 

0209 

PRINT  297,  (FlXtj), J*1,N> 

00002780 

0205 

2000 

CONTINUE 

00002800 

C 

A  CYCLE  EXAMINES  ALL  THE  FACILITIES. 

00002803 

L 

C 

IF  IN  A  CYCLE,  THE  CAPACITY  RULE  RESULTS  IN  SETTING 

00002810 

C 

ADDITIONAL  CX(I,J»  VALUES  EOUAL  TO  2,  THEN  REPEAT 

00002820 

C 

THE  CYCLE.  BUT  IF  FIXIJ»*l  FOR  ALL  J,  THEN  DO  NOT 

00002830 

C 

REPEAT  THE  CYCLE. 

00002835 

C 

0206 

IF  <LQ1.E0.N»  GO  TO  2300 

00002890 

0207 

IF  (LU2.EQ.LR21  GO  TO  2300 

00002895 

0208 

2200 

LR  2*LQ2 

00002860 

L 

0209 

GO  TO  300 

00002870 

C 

*****•••*•***• SOLVE  (L AGRANGI AN  1  RELAXED  PROBLEM************** 

00002680 

C 

FINO  FLBCK1  —  VECTOR  OF  FAC1LIES  FOR  COMPUTING 

00002890 

C 

C 

01,3)  MATRIX  L  LOWER  BOUND.  IT  HAS  VALUE  1  IF  A 

00002900 

c 

FACILITY  IS  USED,  OTHERWISE  IT  HAS  0  VALUE. 

00002910 

0210 

2300 

DO  2500  K* 1 , P 

00002920 

C 

0211 

FLB(KI*0 

00002930 

0212 

2500 

CONTINUE 

000029«»0 

0213 

DO  3000  J*1,N 

00002950 

0219 

IF  (FIX(J). EO.O)  GO  TO  3000 

00002960 

0215 

INO 1 *F1XI ( J ) 

00002970 

L 

0216 

DO  2550  K*1,P 

00002990 

0217 

IF  (EtlNDI.M.EU.O)  GO  TO  2550 

00003000 

0218 

IF  IFLBIM.EO.l)  GO  TO  2550 

00003010 

c 

0219 

FLB(K)*1 

00003020 

0220 

2550 

CONTINUE 

00003030 

0221 

3000 

CONTINUE 

00003060 

l 

0222 

IF  ( I STEP.NE .2 )  GO  TO  3150 

00003070 

0223 

PRINT  3100,  <FLB(K),R*l,P) 

00003080 

0229 

3100 

F0RMAT('0*,MFLB(M,K*1,P>  •,  20I9/16X.20I9) 

00003090 

C 

COMPUTE  COST  MATRIX  C(I*J)  FOR  THE  RELAXED  PROBLtM 

00003100 

c. 

0225 

3150 

DO  3900  J*1,N 

00003110 

0226 

DO  3300  1*1, M 

00003120 

c 

0227 

BSUM*0.0 

00003130 

0228 

DO  3200  X*1,P 

00003190 

0229 

IF  (FLB<M.EQ.ll  GO  TO  3200 

00003150 

e 

0230 

IF  (EII,KI.EQ.0I  GO  TO  3200 

00003160 

02  31 

6SUM*BSUM*(8(M  •  < FLOAT (D( 1 , J,K  I  >/  FLOAT (SIKI1I) 

00003170 

0232 

3200 

CONTINUE 

00003180 

0233 

3250 

cti, ji*a< i,ji*6 sum 

00003190 

i 

0239 

3300 

CONTINUE 

00003200 

0235 

3900 

CONTINUE 

00003210 

0236 

IF  (ISTEP.NE.2)  GO  TO  3995 

00003220 

iv 

0237 

DO  3930  1*1 ,M 

00003230 

0238 

PRINT  3920,  I,  ( C ( I ,  J  )  «  J*1 ,N ) 

00003250 

0239 

3920 

FORMAT  1 /5X ,  •C(I,J)',5X,  •!*•,  13, 2X,  5F15.9, 

00003260 

] 

1  6I/23X,  5F 15.9 ) I 

00003265 

0290 

3930 

CONTINUE 

00003270 

C 

FINO  SUM  OF  MINIMUM  CII,J>  VALUES  OVER  EACH  J, 

00003290 

1 

C 

I.E.,  MINSC*SUM  OF  M1NCU). 

00003300 

C 

IF  F1X(J)*1,  THEN  M INC ( J ) «C ( I , J )  WHERE  CXII,J)*1 

00003310 

35 


T-423 


i 

[ 

S 

i 

! 


FORTRAN 

IV  G  LEVEL 

21  MAIN  DATE  *  80206  14/35207 

0241 

3445 

HINSOQ.O 

00003320 

0242 

DO  3900  J*1,N 

00003340 

0243 

IF  <F1X(J) .EQ.0)  GO  TO  3500 

00003350 

0244 

1N01*FIXI( j» 

00003360 

0245 

MINC1 J)*CUND1,J) 

00003370 

0246 

S0LX( J)*1N0I 

00003380 

0247 

GO  TO  3850 

00003410 

0248 

3500 

LK*0 

00003430 

1 

0249 

1*1 

00003440 

C 

SKIP  CU.JJ  ELEMENT  IF  CXU,J>*2  t  MOVE  TO  NEXT  I 

0000347P 

0250 

3550 

IF  ICX( I,J) .EQ.2J  GO  TO  3700 

00003480 

0251 

IF  Ul-LKj.EQ.l)  GO  TO  3600 

00003485 

0252 

IF  (CU.JJ. GE.MINCtJJJ  GO  TO  3750 

00003490 

0253 

3600 

MlNC(J)*CII,Jl 

00003500 

s 

0254 

IM1N=1 

00003510 

0255 

GO  TO  3750 

00003520 

0256 

3700 

LK=LK*l 

00003530 

U 

02  57 

3750 

1*1*1 

00003590 

0258 

3800 

IF  U.LE.MJ  GO  70  3550 

00003600 

0259 

SOLX  (  J)  *  1 M I  N 

00003610 

1 

0260 

3850 

MINSC*MINSC*M1NC«  J) 

00003620 

0261 

3900 

CONTINUE 

00003630 

0262 

IF  USTEP.NE.2J  GO  TO  3940 

00003640 

w- 

0263 

00  3920  J*1,N 

00003650 

0264 

PRINT  3910  •  J.MINCUI  .SOLXlJ  J 

00003660 

0265 

3910 

FORMAT  1 •0',*J,MINCC J),S0LXCJJ‘,  IS,F15.4,I6J 

00003670 

1 

0266 

3920 

CONTINUE 

00003680 

C 

COMPUTE  FIXED  COST  FC  FOR  LOWER  BOUND 

00003710 

0267 

3940 

FC*0 

00003720 

V- 

0268 

00  4000  K*1,P 

00003730 

0269 

IF  4FL8 (K).EO.O)  GO  TO  4000 

00003740 

0270 

3950 

FC  =  FC*B ( K  J 

00003750 

0271 

4000 

CONTINUE 

00003760 

C 

00003770 

0272 

4050 

L0LB*M1NSC*FC 

00003780 

0273 

IF  (ISTEP.EO.OJ  GO  TO  4150 

00003790 

0274 

PRINT  4120,  M1NSC,  FC,  LOWB 

00003800 

0275 

4120 

FORMAT  CO*,*  MINSC ,  FC,  LOWB  »,  F15.4,  US,  F15.4) 

00003810 

i 

C 

COMPARE  LOWER  BOUND  WITH  BEST  UPPER  BOUNO  STAR 

00003820 

c 

BUBS  WHICH  EQUALS  BU6/ll*EPSJ.  IF  LOWB  IS 

00003830 

c 

GREATER  THAN  OR  EQUAL  TO  BUBS,  THEN  BACKTRACK 

00003840 

0276 

4150 

IF  (LOWB.GE.BUBSJGO  TO  6200 

00003850 

C 

CHECK  IF  CURRENT  SOLUTION  SATISFIES  CAPACITY 

00003880 

C 

CONSTRAINTS 

00003890 

1 

0277 

“<200 

IF  UUNCAP.EO.l)  GO  TO  4420 

00003900 

0278 

4210 

00  4400  K*1,P 

00003910 

0279 

NSUMD-0 

00003920 

f_ 

0280 

00  4300  J«1,N 

00003930 

0281 

1X*S0LXIJI 

00003950 

0282 

NSUMD=NSUM0*0( IX,  j,  kj 

00003960 

0283 

O 

o 

CONTINUE 

00003970 

0284 

IF  USTEP.NE.2J  GO  TO  4320 

00003980 

0285 

PRINT  4310,  K.NSUMO 

00003990 

0286 

4310 

FORMAT  «,0‘,  •K,NSUMD*,2I10J 

00004000 

0287 

432  0 

if(nsumo.le.sikj)  go  to  4400 

00004010 

0288 

GO  TO  5100 

00004020 

- 

0289 

* 

* 

o 

o 

CONTINUE 

00004030 

c 

00004040 

f 
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0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 


0312 

0313 

0314 

0315 

0316 

0317 

0318 

0319 


0320 

0321 


0322 

0323 

0324 

0325 

0326 

0327 


C  ARE  SATISFIED. 

C  UPB*SUM  OF  AlItslt+FlAEO  COST  FCUB  BASED  ON 

C  SOLUTION  VECTOR  SOLXIJt 

C  VECTOR  OF  FACILITIES  FOR  UPPER  BOUND  FUBIKl  HAS 

C  VALUES  1  OR  0  BASED  ON  FACIT1LY  USED  OR  OTHERWISE 

4420  DO  4450  K* 1  * P 
FUBIK.1«0 
4450  COFTINUE 
N$UMA*0 
FCUB*0 

4500  00  4650  J*1,N 
lXsSOLXIJi 
NSUMA*NSUMA+A( IX,  J) 

4550  00  4600  Ke 1 1  P 

IFI E I IX ,K I .EQ.O I  GO  TO  4600 
lFIFUBIM.EQ.il  GO  TO  4600 
FUBIK i*l 
FCU8*FCuB*BIKI 
4600  CONTINUE 
4650  CONTINUE 

IF  (ISTEP.NE.2I  GO  TO  4700 
PRINT  4660,  ( FUB ( K I ,  K— I ,  P 1 
4660  FORMAT  I • 0 • » • I FuB I K I » K* 1 , P )  *,  20I4/16X, 201 4 1 

4700  UPB-NSUMAyFCUB 
4708  IF  IlSTEP.EQ.O)  GO  TO  4750 

PRINT  4710,  NSUMA,  FCUB,  UPB,  BUB,  BUBS 
4710  FORMAT I  ' 0* «  *NSUMA ,  FCUB,  UPS,  BUS,  BUBS  *,2»10,  3F15.4) 

C  COMPARE  UPPER  BOUND  WITH  BEST  UPPER  BOUND 

C  IF  UPB  IS  LESS  THAN  BUB,  SET  IT  aS  BUB  AND 

C  NOTE  THE  SOLUTION 

4750  IF  (UPB.GE.BU6)  GO  TO  5100 
4770  BUB *UPB 

BUB  S=  BUB/  I 1.0+EPS) 

I8N00=N00 
DO  4800  J*1,N 


4800  B50LXI JI*SOLX| J1 
00  4850  X*  1  *  P 
4850  BSOLVIM=FuBIK> 

C  ***«**»*******cOMPARE  LOWB  WITH  8 UBS.  IF  LOWS  IS  GREATER 

C  THAN  OR  EQUAL  TO  BUBS,  THEN  8 ACKTRACK 

4900  IF  ILOWB .GE.BuBSIGO  TO  62U0 

C  «**»*»»*****»*1F  FIXIJl  VALUES  ARE  l  FOR  EACH  J,  THEN  BACKTRACK 

5100  IF  ILQ1.EQ.N)  GO  TO  6200 

C  ***»****»»*»•* APPLY  THE  0OUNOING  RoLE***9****»»»**P9********** 

C  IF  THE  DIFFERENCE  BETWEEN  CII.JI  AND  MINCIJl  IS 

C  GREATER  THAN  THE  DIFFERENCE  BETWEEN  BUBS  AND 

C  LOWB,  THEN  CXI  I , J 1  =  2 

c  *»*,****. *,»,*applY  BRANCHING  RULE  ANO  FIND  BRlt  THE  NEXT 

C  VARIABLE  FOR  LEFT  BRANCHING. 

C  FIND  NMINCIJI,  THE  NtXT  HIGHER  VALUE  THAN  MINCIJl 

C  AND  DIF8RIJI,  THE  DIFFERENCE  BETWEEN  THEN. 


DB0UN0=BU8  S-L0W8 
5200  DO  5250  J«1,N 

NHINC I J I =0. 0 
OIFBR I J 1*0.0 
5250  CONTINUE 

DO  5600  J«1,N 

C  SKIP  TO  NEXT  J  IF  FIXIJl-l 


00004050 

00004060 

00004070 

00004080 

00004090 

00004100 

00004110 

00004120 

00004130 

00004140 

00004150 

00004170 

00004180 

00004190 

00004200 

00004210 

00004220 

00004230 

00004240 

00004250 

00004260 

00004270 

00004280 

00004290 

00004300 

00004310 

00004320 

00004330 

00004340 

00004350 

00004360 

00004370 

00004380 

00004385 

00004390 

00004400 

00004410 

00004420 

00004430 

00004440 

00004450 

00004480 

00004500 

00004510 

00004515 

00004520 

00004525 

00004530 

00004540 

00004550 

00004555 

00004568 

00004570 

00004580 

00004590 

00004600 

00004610 

00004620 
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0326 

IF  (FIXIJ).EQ.II  GO  TO  5600 

00004630 

0329 

LK*0 

00004640 

0330 

1*1 

00004650 

C 

SKIP  CII.JI  IF  CXI  I ,j 1*2  6 

MOVE 

TO  NEXT  1 

00004670 

0331 

5300 

IF  (CXI l.J) .EQ.2I  GO  TO  5350 

00004680 

0332 

IF  ( I  .E U. SOLA ( J I >  GO  TO  5350 

00004690 

0333 

IF  I (Cl I.Jl-MINCI Jl I.GT.DBOUND)  GO  TO 

5330 

00004700 

0334 

IF  (II-LKl.EQ.il  GO  TO  5320 

00004710 

0335 

IF  (Cll.JI .GE.NMINCIJI)  GO  TO  5400 

00004720 

0336 

5320 

NM INC  I J  l*C  1 1  ,J  1 

00004730 

0337 

GO  TO  5400 

00004735 

0336 

5330 

CXI  1 1 j) *2 

00004740 

0339 

NSX=NSX*1 

00004742 

0340 

STX(NSX>*-{ I*100*J|-1000000 

00004745 

0341 

KT2(JI*kT2(J)+1 

00004747 

0342 

IF(KT2(J).LT . ( M— 1 1 1  GO  TO  5350 

00004750 

0343 

IND I  *  SO  LX  I  J  I 

00004752 

0344 

CXI INDI t Jl*l 

00004755 

0345 

NSX*NSX«1 

00004758 

0346 

STXINSXI*  I INDI*100+J>*1000000 

00004760 

0347 

FIXIJ 1*1 

00004762 

0348 

LQ1*LQ1*1 

00004764 

0349 

F IX 1 1  J)  *INO  I 

00004766 

0350 

GO  TO  5600 

00004768 

0351 

5350 

LK*LK*1 

00004770 

0352 

5400 

1*1*1 

00004775 

0353 

IFII.LE.MI  GO  TO  5300 

00004780 

0354 

5500 

D1FBR  I J  l*NM  INC!  J I  «*M  INC  ( J  I 

00004785 

0355 

5600 

CONTINUE 

00004790 

0356 

IF 

1 1 STEP. NE . 2 1  GO  TO  5650 

00004795 

0357 

DO 

5620  1*1, M 

00004820 

0358 

PRINT  290,  I,ICXI1,J),J*1,NI 

00004830 

0359 

5620 

CONTINUE 

00004850 

0360 

PRINT  297,  (FIXIJ), J*1,N) 

00004860 

C 

IF  FIXI J  1*1  FOR  ALL  J,  THEN 

BACKTRACK. 

00004880 

0361 

5650 

IF 

(LQ1.EQ.NI  GO  TO  6200 

00004890 

C 

FINO  MAXDIF,  THE  MAXIMUM  DIFFERENCE  DIFBRIJl 

00004900 

0362 

LF*0 

00004905 

0363 

00 

5800  J*1,N 

00004910 

0364 

IF  (FlXIJl.EQ.il  GO  TO  5690 

00004915 

0365 

IF  I  1 J—LF I . EQ. I  I  GO  TO  5660 

00004920 

0366 

IF  (DIFBRIJl. LT. MAXDIF}  GO  TO  5600 

00004925 

0367 

5660 

MAX0IF=DIFBR|J| 

00004930 

0368 

L  J*  J 

00004935 

0369 

GO  TO  5800 

00004940 

0370 

5690 

LF  *LF*1 

00004943 

0371 

5800 

CONTINUE 

00004946 

0372 

IF 

I ISTEP.N6.2 I  GO  TO  5840 

00004950 

0373 

DO 

5820  J* 1 »N 

00004953 

0374 

IF  (FIXIJI.EQ.lt  GO  TO  5820 

00004956 

0375 

PRINT  5610,  J,  NMINCUI,  MINCIJI,  DIFBRIJl 

0000*960 

0376 

5810 

FORMAT  I *0 • *  •  J .NMI NC  1 J  I » M I NC  C  J ItDlEBRIJI* 

.  15, 

3FI5.4I 

00004963 

0377 

5820 

CONTINUE 

00004966 

C 

••••••»«******bRANCHING  VARIABLE  BR1  CORRESPONDS  TO  MAXDIF***** 

00004970 

0378 

5840 

DO 

5900  J-l.N 

00004980 

0379 

IF  (J.NE.LJ)  GO  TO  5900 

00004990 

0380 

5850 

BR l *SOLX( J  1  *I00*J 

00005000 

0381 

IF 

1 1ST EP. EQ.OI  GO  TO  6020 

00006010 
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0382 

0383 

0386 

0389 


0386 

0387 

0388 

0389 

0390 


0391 

0392 

0393 

0396 

0393 

0396 

0397 

0398 

0399 


0600 

0601 

0602 

0603 

0606 

0605 

0606 

0607 

0608 

0609 

0610 

0611 

0612 

0613 

0616 

0615 

0616 

0617 

0618 

0619 

0620 

0621 

0622 


0623 

0626 

0625 

0626 

0627 

0628 

0629 


PRINT  5880*  BR1 
5880  FORMAT!  •0*«*  BR1SI10I 
GO  TO  6020 
5900  CONTINUE 

C  **************uf>OATE  STX(INS)  AND  NSX*******************.**.** 

C  NSX  REPRESENTS  THE  NUMBER  OF  VARIABLES  IN  STXIlNS) 

6020  NSX*NSX*1 
6060  STX(NSX)*6R1 

IF  I1STEP.NE.2)  GO  TO  6100 
PRINT  6088*  (STX(INS),  INS-1. NSXI 
6088  FORMAT  t  *0*  *  •  STXIINSH,  10110,  1221/,  lOX.lOIlOl) 
c  ******** ****** MOVE  To  The  NEXT  ileft  BRANCH)  nooe  a no  apply 
c  CAPACITY  RULE 

6100  N00-N00*1 

6110  IF  (ET .60.0.0)  GO  TO  6150 
IF ( INSET. EQ.  1 )  GO  TO  6167 
IF  (lNET.EU.l)  GO  TO  6150 
CALL  TIMET  I INT ) 

ELTN*( 1NT-IT0I V26.06E-6 
IF  ltLTN.LT.ET)  GO  TO  6150 
6120  PRINT  6125*  NOD,  ELTN*  BUB*  BUBS*  1BN00 

6125  FORMAT  CO1,  *WAS  AT  N00E**I6*  »  AT  ELAPSED  TIME  **,  F10.6, 

1  •  SECONDS.  S/IX,  *  8UB-SF15.6,  •  BUBS-1  .F15.6, 

2  *  AT  N00E-SI7) 

IBUB-BUB 

IF  I 1BUB.E0. 9999999)  GO  TO  6166 
6130  PRINT  6135*  IB  SOLX I J ) , J-l. N) 

6135  F0RMATI*0S  • SOLUTION  CORRESPONDING  TO  BUB  IS*,  //IX, 

I  MBSOLXU),  J*1,N)S  1018, 3I/18X,  1018)) 

6160  PRINT  6165,  (BSOLY (K ) ,  R«1,P) 

61N5  FORMAT(/lX,'(BSOLY(K),  K*1,P)*,10I8,2I/1BX,  1018)) 

6166  INET-1 
IN1S-ISTEP 
INSET* 1 
ISTEP-2 

GO  TO  6150 

6167  1STEP-1N1S 
INSET-0 

6150  GO  TO  272 

C  ************** ENO  IF  AT  THE  ROOT  NODE***** ***••**•**•*• **••*••• 
6200  IF  INSX.EQ.O)  GO  TO  8100 

6220  IF  I  IABSISTXI NSX) ).GT. 1000000)  GO  TO  6500 
6250  BRO-STXINSX) 

6270  STX(NSX)— BROMOOOOQO 

IF  I1STEP.EQ.0)  GO  TO  6308 
PRINT  6305,  BRO 
6305  FORMAT! *0*,  *BRO  *,110) 

6308  IF  ! 1 STEP.NE. 2 )  GO  TO  6330 

PRINT  6088,  ISTX(lNS),  1NS-1.NSX) 

C  ***<  •««*,, ,,,,MOve  T0  Thfc  NExT  I R IGHT  BRANCH)  NOOE  AND  APPLY 

C  CAPACITY  RULE 

6330  N0D*N00»1 

6610  IF  ItT .EO. 0.0)  GO  TO  6650 
IF  (INSET.tU.ll  GO  TO  6665 
IF  IINET.EU.1)  GO  TO  6650 
CALL  TIMET  I INT I 
ELTNM lNT-lT0»*26.06E-6 
IF  (ELTN.LI.ETI  GO  TO  6650 


00005020 

00005030 

00005060 

00005050 

00005060 

00005070 

00005090 

00005100 

00005150 

00005160 

00005170 

00005220 

00005230 

00005260 

000052*2* 

0000526** 

000052*6* 

00005268* 

00005250* 

00005253* 

00005256* 

00005260* 

00005263* 

00005266* 

00005267* 

00005268* 

00005*70* 

00005273* 

00005276* 

00005280* 

00005290* 

00005292* 

00005296* 

00005296* 

00005298* 

00005300* 

00005302  * 

00005306* 

00005306 

00005308 

00005310 

00005320 

00005330 

000053*0 

00005390 

00005*00 

00005*10 

00005*20 

00005*30 

00005*90 

00005500 

00005510 

00005512* 

00005516* 

00005518* 

00005520* 

00005523* 

00005526* 
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0630 

6620 

PRINT  6125 •  NODt  ELTN,  BUB,  BUBS,  IBNOO 

00005528* 

0631 

IBUBaBUB 

00005530* 

0632 

IF  (IBUB.E0.9999999)  GO  TO  6662 

00005532* 

0633 

6630 

PRINT  6135,  IBS0LXIJ),J*1,N) 

00005533* 

0636 

6660 

PRINT  6165,  (BSOLYU),  *«1,P) 

00005536* 

0635 

6662 

INET*l 

00005538* 

0636 

INIS-1STEP 

00005560* 

0637 

INSET*! 

00005562* 

0638 

i$TEP*2 

00005566* 

0639 

GO  TO  6650 

00005566* 

0660 

6665 

I STEP*1N1S 

00005568* 

0661 

INSET«0 

00005550* 

0662 

6650 

GO  TO  210 

00005552 

0663 

6500 

IF  1  STXINSX).GT. 1000000)  GO  TO  6520 

00005555 

0666 

LX*— ST XI  NSX 1-1000000 

00005560 

0665 

IX*LX/100 

00005570 

0666 

JX*LX— 1X*100 

00005580 

0667 

CX'1X,JX)*0 

00005590 

0668 

KT2 ! JX ) *KT2 1 JX )— 1 

00005595 

0669 

GO  TO  6550 

00005600 

0650 

6520 

LX*  STX ( NSX )— 1 000000 

00005610 

0651 

IX*LX/100 

00005620 

0652 

JX*LX-IX»100 

00005630 

0653 

CXIIX,JX>*0 

00005660 

0656 

FIX) JX)*0 

00005650 

0655 

LOl-LQl-1 

00005660 

0656 

6550 

NSX*NSX— 1 

00005650 

0657 

GO  TO  6200 

00005700 

C 

00005730 

0658 

8100 

IBU8*8UB 

00005760 

0659 

CALL  TIMETIITl) 

00005750* 

0660 

ELT 1*( IT  l-IT0)*26.06E-6 

00005760* 

0661 

PRINT  8105,  ELT1 

00005770* 

0662 

8105 

FORMAT  ('O', ///IX,  'ELAPSED  TIME  IN  SECONDS*',  F15.8) 

00005780* 

0663 

PRINT  8120,  NOD 

00005790 

0666 

8120 

FORMAT  I  'O', 'TOTAL  NUMBER  OF  NODES  EXPLORED  *',I3) 

00005000 

0665 

IF  I IBUB.EU. 9999999)  GO  TO  8350 

00005810 

0666 

8130 

PRINT  8150 

00005820 

0667 

8150 

FORMAT  ('O',  'NOTE)  1.  FOLLOWING  XII ,3)  VARIABLES  SHOW  DESIGN', 

00005830 

1 

l  •  I  TO  WHICH  ACTIVITY  J  IS  ASSIGNED  FOR  J*1  TO  N.', 

00005860 

2 

’  /7X,  '2.  IF  EPSILON  EPS  WAS  ASSIGNED  A  POSITIVE', 

00005850 

3 

1  •  IN0N-2ER0)  VALUE,  THE  SOLUTION  MAY  BE  SUBOPTIMAL. • ./ I 

00005860 

0668 

8180 

PRINT  8200,  IBS0LXI3 ),J*1,N) 

00005870 

0669 

8200 

1 

FORMAT! 'O', T55,  'OPTIMAL  SOLUT ION* ,/lX,T 55, 

00005880 

00005890 

2 

3 1/  ,21X  ,  1018  >  ) 

00005900 

0670 

8220 

PRINT  8250,  IBSOLYIK),  K*1,P) 

00005910 

0671 

8250 

FORMAT  I'O',  • Yl R)  VALUES:',  8X,  1018,  21 / ,21X, 10 18 )  ) 

00005920 

0672 

8280 

PRlIjT  8300, 1BUB 

00005930 

0673 

8300 

FORMAT  ('O',  'OPTIMAL  VALUE  OF  0B3ECT I VE  FUNCTION:',  115///) 

00005960 

0676 

GO  TO  8500 

00005950 

0675 

8350 

PRINT  8600 

00005960 

0676 

8600 

FORMAT  I'O',  •  PROBLEM  DOES  NOT  HAVE  A  FEASIBLE  SOLUTION', 

00005970 

1 

/IX,  •  BECAUSE  THE  CAPACITY  CONSTRAINTS  CANNOT', 

00005980 

2 

/IX,  •  BE  SATISFIED.',/) 

00005990 

0677 

8500 

PRINT  8550 

00006000 

0678 

8550 

FORMAT  1 '0','**6*aN0RMAL  ENO  OF  JOB*****',/) 

00006010 

0679 

8600 

STOP 

00006020 

0680 

ENO 

00006030 

*Relevant  to  the  TIMET  Subroutine 
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SYMBOLIC 

NAME 


DEFINITION 


A(I,J) 

B(K) 

BRO 

BR1 

BSOLX(J) 

BSOLY(K) 

BSUM 

BUB 

BUBS 

C(I,J) 

CX(I,J) 

D(I,J,K) 

DBOUND 

DIFBR(J) 

E(I,K) 

ELT1 


VARIABLE  COST  OF  ACTIVITY  J  USING  DESIGN  I 
FIXED  COST  OF  FACILITY  K 

BRANCHING  VARIABLE  WITH  VALUE  0,  I.E.,  RIGHT-BRANCH 
VARIABLE  X(I,J) 

BRANCHING  VARIABLE  WITH  VALUE  1,  I.E.,  LEFT-BRANCH 
VARIABLE  X(I,J) 

CURRENT  BEST  SOLUTION  CONTAINING  X(I,J)  VARIABLES 
CURRENT  BEST  SOLUTION  CONTAINING  Y (K)  VARIABLES 
B(K)  *  D(I ,  J,  K)/S(K)  SUMMED  OVER  APPLICABLE  K 
BEST  UPPER  BOUND 

BEST  UPPER  BOUND  IF  SOLUTION  MAY  BE  SUBOPTIMAL  WITHIN 
EPS.  BUBS  =  BUB/ (1  +  EPS) 

COST  ELEMENTS  FOR  THE  RELAXED  PROBLEM 

CONTAINS  FIXED  AND  FREE  X(I,J)  VARIABLES.  FIXED 
VARIABLES  ARE  ASSIGNED  A  VALUE  OF  1  OR  2,  AND  FREE 
VARIABLES  ARE  ASSIGNED  A  VALUE  0.  A  VALUE  OF  1  IMPLIES 
THAT  THAT  PARTICULAR  X(1,J)  HAS  A  VALUE  1,  AND  FIX(J)  IS 
FIXED  AS  1.  A  VALUE  OF  2  IMPLIES  THAT  THAT  PARTICULAR 
VARIABLE  SHOULD  NOT  BE  CONSIDERED  FOR  CURRENT  COMPUTATIONS 
CX(I,J)  IS  UPDATED  BY  THE  CAPACITY  RULE,  THE  BOUNDING  RULE, 
AND  THE  RULES  FOR  BRANCHING  AND  BACKTRACKING 

CAPACITY  REQUIRED  AT  FACILITY  K  FOR  ACTIVITY  J  WHEN 
ACTIVITY  J  USES  DESIGN  I 

DIFFERENCE  BETWEEN  BOUNDS  BUBS  AND  LOWB 

DIFFERENCE  IN  C(I,J)  VALUES  FOR  SELECTING  BRANCHING 
VARIABLE 

HAS  VALUE  1  IF  DESIGN  I  USES  FACILITY  K,  OTHERWISE 
THE  VALUE  IS  0 

ELAPSED  TIME  IN  SECONDS  TO  EXECUTE  THE  PROGRAM  OBTAINED 
FROM  THE  TIMET  SUBROUTINE 


ELTN 


ELAPSED  TIME  IN  SECONDS  TO  COMPARE  WITH  ET  IF  SPECIFIED 


EPS 

ET 

FC 

FCUB 

FIX(J) 

FIXI(J) 

FLB(K) 

FUB(K) 

I 

IBALD 

IBNOD 

IBUB 

ICAPR 

IDIFD 

IINPT 

IMIN 

INET 

INIS 

INS 

INSET 

INT 


EPSILON  VALUE  IF  A  SUBOPTIMAL  SOLUTION  ACCEPTABLE,  E.G, 
EPS  OF  0.002  IMPLIES  THAT  THE  SOLUTION  MAY  BE  SUBOPTIMAL, 
BUT  IS  GUARANTEED  TO  HAVE  A  VALUE  WITHIN  0.2%  OF  THE 
OPTIMAL  VALUE.  EPS  IS  0  IF  AN  OPTIMAL  SOLUTION  IS 
DESIRED 

ELAPSED  TIME,  IF  SPECIFIED,  AT  WHICH  THE  NODE  AND  BOUNDS 
INFORMATION  IS  PRINTED 

FIXED  COST  FOR  COMPUTING  LOWER  BOUND 

FIXED  COST  FOR  COMPUTING  UPPER  BOUND 

VECTOR  OF  FIXED  COLUMNS.  A  COLUMN  IS  FIXED  IF  IT  HAS 
A  FIXED  X(I,J)  VARIABLE  WITH  VALUE  1 

VECTOR  SPECIFYING  INDEX  I  CORRESPONDING  TO  FIX(J)  OF 
VALUE  1 

VECTOR  OF  FACILITIES  FOR  COMPUTING  LOWER  BOUND 
VECTOR  OF  FACILITIES  FOR  COMPUTING  UPPER  BOUND 
INDEX  FOR  DESIGNS 

BALANCE  AVAILABLE  CAPACITY  FOR  A  GIVEN  FACILITY 
NODE  CORRESPONDING  TO  THE  BEST  UPPER  BOUND 
BEST  UPFER  BOUND 
OPTION  TO  USE  CAPACITY  RULE 

DIFFERENCE  BETWEEN  A  D(I,J,K)  VALUE  AND  MIND(J) 

OPTION  TO  LIST  INPUT  DATA 

MINIMUM  OVER  INDEX  I  FOR  A  GIVEN  COLUMN  J 

INDICATOR  FOR  ET 

HAS  THE  SAME  VALUE  AS  I STEP 

INDEX  FOR  STACK  OF  VARIABLES  STX 

INDEX  FOR  DISPLAYING  DETAILED  STEPS  IF  OPTION  ET  IS 
SPECIFIED 

INTERMEDIATE  CALL  TO  TIMET  SUBROUTINE 


T-423 


ISTEP 

ITO 

IT1 

I  UN  CAP 

IX 

J 

JX 

K 

KT2 (J) 

LK 

LOWB 

LQ1 

LQ2 

LR2 

P 

S(K) 

SOLX(J) 

STX(INS) 

TIMET 


OPTION  TO  LIST  SUMMARY  OF  INTERMEDIATE  STEPS,  DETAILED 
STEPS,  OR  TO  SKIP  THE  LISTING 

INITIAL  CALL  TO  TIMET  SUBROUTINE 

FINAL  CALL  TO  TIMET  SUBROUTINE 

OPTION  TO  INDICATE  IF  THE  PROBLEM  BEING  SOLVED  IS 
CAPACITATED  OR  UNCAPACITATED 

INDEX  I  OF  VARIABLE  X(I,J) 

INDEX  FOR  ACTIVITIES 

INDEX  J  OF  VARIABLE  X(I,J) 

INDEX  FOR  FACILITIES 

A  COUNTER  FOR  CX(I,J)  VALUES  SET  EQUAL  TO  2  FOR 
COLUMN  J 

A  COUNTER 

LOWER  BOUND 

A  COUNTER  FOR  CX(I,J)  VALUES  SET  EQUAL  TO  1 

A  COUNTER  FOR  CX(I,J)  VALUES  SET  EQUAL  TO  2  IN  A 
CYCLE  IN  WHICH  ALL  THE  FACILITIES  ARE  EXAMINED 

A  COUNTER  SET  EQUAL  TO  LQ2  AT  THE  END  OF  A  CYCLE  WHEN 
APPLYING  THE  CAPACITY  RULE 

NUMBER  OF  FACILITIES 

TOTAL  CAPACITY  OF  FACILITY  K 

SOLUTION  SHOWING  INDEX  1  OF  X(I,J)  VARIABLES 
CORRESPONDING  TO  J=l,  2 . N 

STACK  OF  FIXED  X(I,J)  VARIABLES  FOR  BRANCHING  AND  BACK¬ 
TRACKING.  AN  X(I,J)  RECORDED  IN  CX(I,J)  AS  1  DUE 
TO  THE  BRANCHING  RULE  IS  RECORDED  IN  STX(INS)  AS 
X  *  100  +  J  .  AN  X(I,J)  RECORDED  IN  CX(I,J)  AS 
1  DUE  TO  THE  CAPACITY  RULE  OR  THE  BOUNDING  RULE  IS 
RECORDED  IN  STX(INS)  AS  (X  *  100  +  J)  +  1,000,000. 

AN  X(I,J)  RECORDED  IN  CX(I,J)  AS  2  IS  RECORDED 
IN  STX(INS)  AS  -(X  *  100  +  J)  -  1,000,000 

A  SUBROUTINE  FROM  THE  PL1  LIBRARY  TO  RECORD  THE 
EXECUTION  TIME 
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UPB  UPPER  BOUND 

ZIPCAP  NAME  OF  THE  PROGRAM,  AN  ACRONYM  FOR  ZERO-ONE  INTEGER 

PROGRAM  TO  SOLVE  MULTIACTIVITY  MULTIFACILITY  CAPACITY- 
CONSTRAINED  ASSIGNMENT  PROBLEMS 
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APPENDIX  C 

ANNOTATED  OUTPUT  FOR  TEST  PROBLEM  1 


0*0  *13  0*0  *?43  0*4V3W«I  1*431?  I  1**401  1«14Nll  l  0313313?  ?N01140 


001  ?i.l  SC  *«1 


001 


s*m  *o  Auburn 


nl*SC>  K ,  LOWS  33328.3633  1730  37278.3633 


orim»L  v»lui  07  osjtcTivt  function*  *oi76 


APPENDIX  D 

ANNOTATED  OUTPUT  FOR  TEST  PROBLEM  2 


50  - 


0*1  IONS  SELECTED  I  I1NPT-1  1CAPM1  IS11P-0  1UNCAP-0  tPS*  0.00200  fcT*  0.0 
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CAPACITY  USAGE  MA1A1X  Oll.J.At 
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APPENDIX  E 

ANNOTATED  OUTPUT  FOR  TEST  PROBLEM  3 
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OPT! QMS  SELECTED  >  IIKTT-I  ICAPt.Q  1STEP-0  I UNCAP* 1  EPS*  0.0  ET-  10.000 


EMEO  COST  VECIWI  81  ft  I 


ctfitin  LiMi?  mtot  sui 
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